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ABSTRACT 


We report on a Spitzer /IRB (mid-infrared), /?ATi?/PCA-|-HEXTE (X- 
ray), and Ryle (radio) simultaneous multi- wavelength study of the microquasar 
Cygnus X-1, which aimed at an investigation of the origin of its mid- infrared 
emission. Compact jets were present in two out of three observations, and we 
show that they strongly contribute to the mid-infrared continuum. During the 
first observation, we detect the spectral break — where the transition from the 
optically thick to the optically thin regime takes place — at about 2.9 x 10^® Hz. 
We then show that the jet’s optically thin synchrotron emission accounts for 
the Cygnus X-l’s emission beyond 400 keV, although it cannot alone explain 
its 3 — 200 keV continuum. A compact jet was also present during the second 
observation, but we do not detect the break, since it has likely shifted to higher 
frequencies. In contrast, the compact jet was absent during the last observation, 
and we show that the 5 — 30 pm mid-infrared continuum of Cygnus X-1 stems 
from the blue supergiant companion star HD 226868. Indeed, the emission can 
then be understood as the combination of the photospheric Raleigh-Jeans tail and 
the bremsstrahlung from the expanding stellar wind. Moreover, the stellar wind 
is found to be clumpy, with a filling factor « 0.09 — 0.10. Its bremsstrahlung 
emission is likely anti-correlated to the soft X-ray emission, suggesting an anti- 
correlation between the mass-loss and mass-accretion rates. Nevertheless, we do 
not detect any mid-infared spectroscopic evidence of interaction between the jets 
and the Cygnus X-l’s environment and/or companion star’s stellar wind. 

Subject headings: binaries: close — X-rays: binaries — Infrared: stars — dust, 
extinction — Stars: individual: Cygnus X-1 — ISM: jets and outflows 



1. Introduction 


Accretion— ejection phenomena are responsible for the strong radio/infrared/X-rays 
correlations observed in microquasars (see e.g. Mirabel et al. 1998; Corbel & Fender 2002; 
Gallo et al. 2003; Coriat et al. 2009). Indeed, some of the accreted material is returned 
to the interstellar medium (ISM) through collimated relativistic ejections, whose power 
law-like synchrotron radiation {F^ oc is thought to extend from the radio to the X-ray 
domains. The best known are the so-called compact jets, only detected in the hard state 
(HS), whose spectral continuum is well-modeled by a flat or weakly-inverted power law in 
the optically thick regime (0 < a < 0.2) and decreasing in the optically thin one, with 
a ranging between —0.4 and —1 depending on the electron energy distribution (for a 
—0.6 canonical value, Blandford & Koenigl 1979; Falcke & Biermann 1995). To date, the 
spectral break where this transition takes place has only been detected photometrically in 
two sources, GX 339—4 and 4U 0614-H091 (Corbel <fe Fender 2002; Migliari et al. 2006, 
2010). Yet the knowledge of this cut-off frequency is crucial for a full understanding of the 
properties of the compact jets, as it is closely related to the original physical conditions, 
such as the magnetic field, the base radius of the jet, and the total energy of the electron 
population. 

Among microquasars, Cygnus X-1 (Bowyer et al. 1965) is one of the few classified as a 
high mass X-ray binary (HAIXB). It is composed of an 8—16 Mq black hole (Shaposhnikov 
& Titarchuk 2007; Caballero-Nieves et al. 2009) bound to the 17—31 Mq 09.7 ab supergiant 
HD 226868 (Walborn 1973; lorio 2008; Caballero-Nieves et al. 2009). It has a quasi-circular 
5.6 d orbit (Brocksopp et al. 1999a, and reference therein), with an associated inclination 
in the range of 23° — 43° (Gies et al. 2003; Caballero-Nieves et al. 2009). The distance to 
the source is still a matter of debate, but several studies have found values between 1.5 and 
2.5 kpc (Ziolkowski 2005; Caballero-Nieves et al. 2009, and reference therein). Recently, 



Xiang et al. (2011) found 1.80 kpc with a 7% accuracy from X-ray dust scattering halo 
studies. Moreover, Reid et al. (2011) derived 1.86 kpc with a 6% accuracy, from radio 
parallax measurements obtained with the very large baseline array (VLBA) at 8.4 GHz 
between 2009 Jan 23 and 2010 Jan 25. Both results are consistent with 2 kpc within the 
uncertainties, and we adopt this value in this study. 

Cygnus X-1 spends two-thirds of its time in the HS (Wilms et al. 2006). Indeed, 
although the companion is close to filling its Roche lobe, the transfer of material takes 
place through the strong stellar winds that are focused towards the black hole (see Gies & 
Bolton 1986; Gies et al. 2003; Hanke et al. 2009). This lower accretion rate transfer is less 
effective than for Roche lobe overflow microquasars, thereby most of the time preventing 
the transition to a disk-dominated soft state (SS, see e.g. Pottschmidt et al. 2000, 2003). 
Instead, the source often exhibits a so-called failed-state transition, characterized by a 
continuous increase of the X-ray flux up to an intermediate state (IS), followed by a strong 
decrease and a return to the HS, with possible launches of discrete ejections. Finally, its 
emission appears to be in phase with the orbital period in most of the spectral domains in 
which the source has been detected, (Pooley et al. 1999; Brocksopp et al. 1999a; Lachowicz 
et al. 2006). This includes the radio activity, which is similar to that of other microquasars. 
In particular, the compact jets, characterized by a very flat spectrum (a < 0.1) and 
a 10—25 mJy flux density level are detected in the HS (Fender et al. 2000; Stirling et al. 2001). 


In this paper, we report on a mid-infrared (MIR), radio, and X-ray multi- wavelength 
study of Cygnus X-1, with a focus on MIR spectroscopy. Our goals are to (1) assess the 
origin of the MIR emission, (2) compare the possible variation of the MIR continuum to 
that of the X-ray and/or radio emission, and (3) detect any spectroscopic feature arising 
from an interaction between the jets, the stellar wind, and/or the interstellar medium 
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(ISM). The observations and the data reduction are presented in Sect. 2, while Sect. 3 
is devoted to the X-ray and radio study. In Sect. 4 and Sect. 5, we analyze the broad 
band MIR to radio SEDs of the source and the spectroscopic contents of the MIR spectra, 
respectively. Finally, we discuss the outcomes in Sect. 6. A companion paper, Lee et al. (in 
prep), will focus on the dust content of Cygnus X-1. 


2. Observations and data reduction 

Cygnus X-1 was observed on UT 2005 May 23, June 7, and July 2 (P.I. S. Heinz) with 
the InfraRed Spectrograph (IRS; Houck et al. 2004) on-board the Spitzer Space Telescope 
{Spitzer: Werner et al. 2004). In the X-ray domain, the source was observed simultaneously 
with the Proportional Counter Array (PCA; Jahoda et al. 2005) and the High Energy X-ray 
Timing Experiment (HEXTE; Rothschild et al. 1998) instruments (P.I. J. Wilms), and 
we made use of the 1.2—12 keV public data provided by the All-Sky Monitor instrument 
(ASM; Levine et al. 1996), all mounted on the Rossi X-ray Timing Explorer (RXTE). 
Finally, Cygnus X-1 was regularly monitored at 15 GHz with the Ryle radio telescope (P.I. 
G. G. Pooley) and we used the 20 s and 5 min light curves that are simultaneous with the 
IRS spectra, for the analysis presented in this paper (see Table 1). 


2.1. IRS observations 

We used the SL2 (5.20 — 7.70 /^m), SLl (7.40 — 14.50 /^m), LL2 (14.00 — 21.30 ^m), and 
LLl (19.50 — 38.00 pm) modules at two nodding positions for sky subtraction. The total 
exposure times were each time set to 60 s in SLl and SL2, 120 s in LL2, and 2x120 s in 
LLl. Basic Calibration Data (BCD) were reduced following the standard procedure given 



in the IRS Data Handbook^. The basic steps include bad pixel correction with IRSCLEAN 
vl.9, sky subtraction, as well as extraction and calibration (wavelength and flux) of the 
spectra — with the Spectroscopic Modeling Analysis and Reduction Tool software SMART 
v8.1.2 — for each nodding position. Note that in LL2 and LLl, where the sky appeared 
weakly non-uniform between both nods, the subtraction residuals were fitted by a first order 
polynomial. Spectra were then nod-averaged to improve the signal-to-noise ratio (SNR). 

The flux calibration appeared consistent between the different modules at Obs. 1 and 
Obs. 2, with differences less than 5%, and this allowed us to combine all the spectra by 
scaling them to SLl, which has the highest flux level with respect to SL2 and LL2. However, 
at Obs. 3, SL2 had a flux level about 14% lower than that of SLl, while no such problem 
was detected with LL2 nor LLl. We were not able to definitively assess the reason for such 
a discrepancy, but a misalignement of the SL2 slit is the most reasonable explanation. We 
therefore scaled the SL2 flux level to that of SLl. 


2.2. RXTE observations 

In the reduction of the data from the two pointed instruments PCA and HEXTE, we 
follow the procedure used in previous studies of Cygnus X-1 (e.g. Wilms et al. 2006). Data 
were extracted using HEASOFT 6.9 for all times at least 10 minutes after passages of the 
South Atlantic Anomaly. During the observations, two (first two RXTE observations) and 
three (third RXTE pointing) proportional counter units were in operation. 


^http://ssc.spitzer.caltech.edu/irs/irsinstrumenthandbook/IRS-Instrument-Handbook.pdf 
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2.3. MIR data dereddening 

Our Cygnus X-1 observations show a strong absorption feature around 9.7 /im. At 
the distance to the source (« 2 kpc), we expect a considerable column of interstellar 
material, so this feature is most likely from silicate dust in the ISM along the line-of-sight. 
Several consistent measurements of the visible extinction are found in the literature (see 
Caballero-Nieves et al. 2009, and reference therein). In the companion paper (Lee et al., 
in prep), we derive an absolute visual extinction Ay « 2.95 from the opacity of the ISM 
silicate absorption feature T9.7, defined by Ay = 18.5 x tq.t (Draine 2003), consistent with 
that of Wu et al. (1982), Ay = 2.95 ± 0.22 with Ry = 3.1. 

We therefore use this value to deredden the spectra with the extinction laws given in 
Chiar & Tielens (2006), one of the most recent works covering the entire IRS spectral range. 
In their paper, the authors derived the Aa/Ak ratio rather than the usual Ax/Ay one. To 
express the extinction ratio in the standard way, we assign to Ak the value derived from the 
law given in Fitzpatrick (1999) for the diffuse ISM {Ry = 3.1), which is Ak = 0.111 x Ay- 
Moreover, these authors find that beyond 8 /um, where silicate absorption dominates, the 
MIR extinction is different for the diffuse ISM and the Galactic center. They consequently 
propose two extinction laws, and we explored both options for dereddening our spectra 
of Cygnus X-1. We find that the Galactic center case clearly overestimates the silicate 
absorption depth, and we therefore adopt the corrections based on the extinction due to the 
diffuse ISM. We note that a single extinction correction (i.e., Ay = 2.95 ± 0.22) produces 
equally good results for all three epochs. This gives us further confidence that the observed 
silicate absorption features are caused by dust in the diffuse ISM as opposed to material 
directly associated with the Cj^gnus X-1 sj^'stem. 

The resulting dereddened spectra displayed in Fig. 1 show' unambiguously that the 
MIR emission of Cygnus X-1 is variable. Indeed, the MIR continuum was brighter during 



Obs. 1 and Obs. 2, the difference with respect to Obs. 3 increasing with wavelength, from 
no difference at 6 fj,m to an excess of 15—20 mJy beyond 10 pm. 


3. X-ray/radio behavior 

Since dust heating as a contribution to the MIR emission can be ruled out due to 
the lack of observed characteristic dust spectral features (Lee et a., in prep), we consider 
four main alternative physical processes: (1) the stellar photosphere, whose emission is 
constant, (2) bremsstrahlung from the expanding stellar wind, (3) synchrotron emission 
from jets, (4) bremsstrahlung from the accretion disk wind. We use the simultaneous radio 
and X-ray data to assess the extent to which the last two components are present during 
the observations (see Sect. 4 for details). 


3.1. RXTE spectra 

The RXTE spectra are fitted between 3 and 200 keV in XSPEC v. 12.5.1. We tested 
several continuum models such as the standard a— disk (Mitsuda et al. 1984) plus a 
Comptonization component (Titarchuk 1994) modified by photoelectric absorption, but we 
find that all the spectra are best described by the phenomenological model combining a 
broken power law with a high energy cut-off and a gaussian to account for the iron emission, 
both modified by photoelectric absorption. These results are consistent with previous 
studies (see e.g. Wilms et al. 2006). The best-fit parameters for all the observations are 
listed in Table 2 and the fits are displayed in Fig. 2. 

Based on the classification of Cygnus X-Ts spectral states as given by Wilms et al. 
(2006), these fit results show unambiguously that the source is in the HS during Obs. 1 
and Obs. 3, wdth Fi photon indices of about 1.83 and 1.75, respectively. The 3 — 200 keV 



unabsorbed flux is nevertheless about 50% larger during Obs. 1 than during Obs. 3 which, 
considering the X-ray/radio correlation observed in microquasars in the HS, suggests a 
stronger contribution of compact jets to the overall emission of Cygnus X-1 during Obs. 1. 
In contrast, the spectrum of the source during Obs. 2 is softer, with Fi « 2.08. Such a 
photon index is consistent with an IS. But similarly, the high 3 — 200 keV unabsorbed flux 
suggests a stronger contribution of compact jets than during Obs. 3. 


3.2. Light curves 

In Fig. 3, we display the Cygnus X-1 radio (15 GHz, 5 min time resolution) and ASM 
(1.2“12 keV, day-averaged) light curves, as well as the ASM HR2 (C/A) hardness ratio, 
between MJD 53500 and MJD 53580, bracketing our observations. Clearly, it is seen that 
the count rate is lower and the hardness ratio higher during Obs. 3, i.e. when the MIR 
continuum appears to be at its lowest level (see Fig. 1, blue spectrum). Consistent with the 
spectral fitting, the ASM light curve and HR2 hardness ratio show that the source during 
Obs. 2 is in a softer state and exhibits a slightly higher flux than during Obs. 1 (see also 
Fig. 4). This is in good agreement with the light curve trend shown in Fig. 3 between 
MJD 53510 and MJD 53530, which shows a sharp decrease in the hardness ratio and a 
slow increase in the count rate. This probably points to a higher accretion disk activity 
becoming unsteady around MJD 53520. Then, the count rate drops between MJD 53530 
and MJD 53550 before stabilizing, while the hardness ratio follow's the opposite path. 
This behavior is characteristics of a failed-state transition as characterized by Wilms et al. 
(2006). 

Moreover, the Ryle light curves during Obs 1. and Obs. 2 exhibit an average radio flux 
between 12 and 25 mJy, while the radio flux is on average lower and less steady during 
Obs. 3. This behavior is shown by the 20 s time resolution Ryle light curves, displayed in 
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Fig. 5, which cover the exact time range of each IRS observation. It is clear that Obs. 1 
and Obs. 2 are consistent with compact jets. In contrast the radio activity during Obs. 3 is 
very unstable, covering a range from non-detection to 48 mJy. This behavior is not by itself 
inconsistent with the presence of compact jets, and possible explanations of this flickering 
could be attributed to some instabilities internal to the jet or some peculiar accretion flow 
properties. 

Finally, it is worth mentioning that there have been indications, from the timing 
behavior, of two different HS regimes (Pottschmidt et al. 2003; Axelsson et al. 2008). One 
of them is characterized by a very high hardness level, a low X-ray flux, and characteristic 
frequencies shifted to the lower edge in the power density spectrum (PDS). This is clearly 
the case of the HS detected during Obs. 3 (see Grinberg et al. in prep, for the full PDS), 
which is actually the beginning of the hardest X-ray activity ever observed for Cygnus X-1 
(Nowak et al. 2011). In contrast, the HS properties during Obs. 1 are consistent with what 
is expected in the canonical picture. 


4. The MIR continuum fitting 

We use spectral fitting to understand the origin of the continuum variation between 
the three IRS observations. Between nod 1 and nod 2, we on average found 4% and 5% 
flux differences in SL and LL modules, respectively. These quantities were quadratically 
added as systematic errors to take into account uncertainties due to flux calibration. Since 
Cygnus X-1 was never in the SS during our observations, we rule out bremsstrahlung from 
the accretion disk wind, and we onty consider thermal emission from the companion star’s 
photosphere, bremsstrahlung from the expanding stellar wind, and/or synchrotron from the 
compact jets. Moreover a bump is present in the three spectra between about 17 and 22 pm. 
The origin of this bump is a matter of debate, and it could be due to the dereddening 
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process or to irradiation of dust around the BH (S. Markoff, private communication), but 
more likely to the Spitzer /IRS> LLl 24 fim deficit^. 


4.1. Obs. 3: the stellar continuum 

Considering the very unstable radio activity of Cygnus X-1 during Obs. 3, a jet 
contribution to the MIR spectrum is very unlikely. We therefore only consider the emission 
from the star. To model the photosphere, we simply use a black body. Cygnus X-l’s 
companion star, HD 226868, is an 09.7ab supergiant and its unabsorbed emission peaks 
in the ultraviolet (UV). The MIR alone is therefore not sufficient to constrain the star’s 
temperature, so we fix it to the most recent derived value, i.e. 28000 K as given in 
Caballero-Nieves et al. (2009), based on theoretical fitting to the ultraviolet spectrum. In 
their seminal papers, Wright &: Barlow (1975) and Panagia & Felli (1975) showed that the 
bremsstrahlung from a homogenous spherically expanding wind would be responsible for a 
MIR and radio excess oc ® in supergiant stars. We therefore add a power law to take the 
wind emission into account. This simple model is summarized by Eq. 1: 

CM =(|)'xB(., 28000) + .-V,(^)- (1) 

where i?* and D* are the companion star’s radius and distance, B is the black body 
emission, and Ng and ag are the bremsstrahlung flux density at 15 GHz and spectral index. 

The best fit is displayed in Fig. 6 and the best-fit parameters are listed in Table 3. 
The ratio is consistent with the one expected from the Cygnus X-l’s companion 

star, since the derived stellar radius for a typical 2 kpc distance is about 19.9 i?©, the 
typical value for a typical late 0 supergiant. Moreover, the power law spectral index is 


"http: II ssc.spitzer.caltech.edu/irs/features/^3-LLl_24uT)eficit 
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characteristics of bremsstrahlung from an expanding stellar wind, with at ~ 0.55, consistent 
with the 0.6 theoretical value within the 90% uncertainties. The MIR continuum of 
Cygnus X-1 during Obs. 3 can therefore be described by that of HD 226868. This alone 
does not rule out the possible presence of compact jets. Nevertheless, both the average 
X-ray and radio fluxes are at least 50% lower than during Obs. 1 and 2. So, even if they are 
present, either the compact jets are too faint to be detected in the MIR, or they break to 
the optically thin regime beyond 30 fim. Indeed, in this regime, the IR flux is expected to 
vary with respect to the radio flux as Fm « Fradio^^^^'^ (Heinz 2004). Considering the low 
radio flux level, this very probably excludes a significant contribution of a compact jet to 
the MIR spectrum of Cygnus X-1 during Obs. 3. 


4.2. Obs. 1 and Obs. 2: MIR contribution from the compact jets 

The radio behavior does exhibit compact jets that may be responsible for the observed 
continuum increase. We consider two cases to model their contribution, (1) the spectral 
break does not occur in the IRS spectral range so the optically thick synchrotron emission 
can be described by a simple power law, or (2) the spectral break is located in the IRS 
spectral range and the emission of the compact jets can be modeled by a broken power law. 

To better constrain the compact jets emission parameters, we include the average 
15 GHz flux densities of Cygnus X-1 during the IRS/LL observations. We then fit these 
MIR/radio spectral energy distributions (SEDs) by fixing the black body parameters and 
bremsstrahlung spectral index to those of Obs. 3. 


4-2.1. Case 1: simple power law 


The overall equation describing the Cygnus X-1 emission is: 
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XB(., 28000) + ,V.(^ 



15 GHz 


(2) 


where iVj and aj are the synchrotron flux density at 15 GHz and spectral index, respectively. 


Here, Ng was left free to mimic a possible change in the bremsstrahlung intensity due 
to a variation of the star’s mass loss rate. The best-fit parameters are listed in Table 4, 
and the corresponding fitted SEDs displayed in Fig. 7. Both fits are consistent with a 
contribution of compact jets to the MIR emission of Cygnus X-1 during Obs. 1 and Obs. 2. 
Indeed, the derived spectral indices are 0.01 ± 0.02 and 0.05 ± 0.01, respectively, consistent 
with typical values. Moreover, the bremsstrahlung flux density at 15 GHz during Obs. 1 is 
found to be similar to that during Obs. 3 (the no-jet case, see Sect. 4.1), 0.66 ± 0.06 mJy 
compared to 0.74io;J° mJy. In contrast, Ng during Obs. 2 is lower (0.441°;“ mJy). 


4-2.2. Case 2: broken power law 
Here, the overall equation describing the Cygnus X-1 emission is: 


F{u) = 



xB{n, 28000) +iVff 


15 GHz 




N, 




J V 15 GHz; 


AT ( Vc / 1/ \«j2 

J Vl5 GHz/ V15 GHz/ 


li V < Vc 
ii V > Vc 


( 3 ) 

where A), aji, a^, and Uc are the jets’ synchrotron flux density at 15 GHz, optically thick 
and thin spectral indices, and cut-off frequency, respectively. 


When wc left free to vary all the parameters of the broken power law^ and the 
bremsstrahlung normalization Ng at the same time; however this results in poorly 
constrained fits due to several degeneracies between the various fit parameters. We therefore 
proceeded wuth three scenarios, each time fixing one parameter only: (a) a ^2 is fixed to —0.6, 
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i.e. the theoretical spectral index expected from the optically thin synchrotron emission of 
the compact jets, (b) Ng is fixed to the values found in case 1, (c) Ng is fixed to the value 
found during Obs. 3, i.e. we assume that the bremsstrahlung intensity is the same for all 
the observations. The best-fit parameters are listed in Tables 5, 6, and 7, and the SEDs are 
displayed in Fig. 8, respectively. 


4-2.3. Interpretation 

A broken power law gives better results to describe the emission of the compact jets 
during Obs. 1, with reduced ^ values of 0.81, 0.85, and 0.78 for scenarios (a), (b), and 
(c), respectively, versus 0.96 for a simple power law. The low reduced ^ are largely 
due to artificially forced high IRS systematic errors that were quadratically added to the 
uncertainties. Nevertheless, A ^ > 34 for only one DOF less, and the three fits give very 
similar results. The spectral break is clustered in the MIR range (2.70 — 2.94) x 10^^ Hz, 
i.e. about 10.06 — 11.11 /rm, and the optically thin spectral index aj 2 for (b) and (c) is 
always consistent with the expected value (between —1 and —0.4) within the uncertainties. 
Moreover, for (a), the bremsstrahlung intensity is found similar to that during Obs. 3, a 
hint that the mass-loss rate is the same during both observations, as expected in the HS. 
So, based on the aforementioned reasons, we conclude that the consistency of all the fits 
strengthens the need for a broken power law to describe the MIR continuum of Cygnus X-1 
during Obs. 1. 

The interpretation of the results for Obs. 2 is more complex, even though the broken 
power law' gives lower reduced ' (0.69 against 0.72). First, A ^ w 9.7, which is far less 
significant than the values found for Obs. 1. Moreover, the fits with the broken powder 
law' are not physically relevant. Indeed, the fit for scenario (b) gives a positive photon 
index after the spectral break {a w 1.21), w'hich is clearly inconsistent with a compact 



jet. The two other fits are consistent with compact jets, but give cut-off frequencies that 
are about twice smaller than during Obs. 1. This is counter-intuitive since the opposite is 
expected. Indeed, the source during Obs. 2 is not in the HS but in a failed-state transition, 
i.e. an IS with a softer and brighter spectrum. In models of optically thick jet emission, 
the break frequency is coupled to the jet luminosity: if the luminosity increases, the break 
frequency also increases. An increase in the optically thick flux requires an increase in the 
jet energy density and/or the cross section of the jet, which always increases the optical 
depth, thus increasing the break frequency. The fact that our best fits for Obs. 2 indicate 
the opposite would imply a change in the jet geometry, e.g., a significant increase in the 
distance of the base of the jet along the jet axis or an increase in the opening angle of the 
jet (such that the footprint of the base of the jet shrinks, decreasing its optical depth). 
Given the level of uncertainty in our data (demonstrated by the fact that our models are 
under-constrained), such conclusions are still rather speculative at this point. Moreover, 
several studies of Cygnus X-1 showed the existence of an anti-correlation between the 
mass-loss and the accretion rates (Gies et al. 2003; Boroson & Dil Vrtilek 2010), maybe 
due to photoionization of the wind by the soft X-ray emission. During the IS and SS, the 
bremsstrahlung contribution should therefore be lower than in the HS. This is actually 
exactly the result of case 1, where we consider a simple power law instead of a broken one. 
The bremsstrahlung intensity is found to be about 40% lower than that during Obs. 3, 
while the cut-off frequency is assumed to occur below .5.48 //m. We believe that it better 
describes what occurs during Obs. 2. 


5. MIR spectral lines analysis 

In the companion paper (Lee et al., in prep), we show that the presence of warm dust 
around Cygnus X-1 is very unlikely and that the strong silicate absorption feature is due 
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to the ISM. Here, we study the spectral line variability of the MlR spectra between the 
three different epochs. Moreover, because the spectral features are remarkably similar for 
the three spectra in the SL modules, we only focus on the longer wavelengths, from 14.5 to 
32 yum. Finally, due to potential rapid variations with respect to the source’s X-ray or radio 
activity, we separately study nod 1 and nod 2, in contrast to using nod-averaged spectra 
for the continuum fitting. All the non-dereddened spectra are displayed in Fig. 9 and the 
features listed in Table 8. Note that we only list the features detected at 5 ct or more above 
the continuum (assessed in 2 /rm windows using a second-order polynomial). 

As seen in Fig. 9 and Table 8, the MIR spectra of Cygnus X-1 are dominated by H I 
and He II emission lines that likely arise from the stellar wind of HD 226868. Overall, more 
features are detected during Obs. 1 and Obs. 2, probably because of a higher soft X-ray 
emission flux level that photoionizes the wind. Moreover, in both cases, the spectroscopic 
content appears to vary from nod to nod, with more lines in nod 1 than in nod 2. This 
is likely due to slight background variations between both positions, although dramatic 
changes of the photoionization state cannot be excluded. 

During Obs. 3, the same features are detected at both positions with the exception of 
two absorption lines at 15.60 ± 0.01 and 25.90 ± 0.02 ^m in the nod 1 spectrum, which 
correspond to [Ne III] and He II. A forbidden absorption feature of Ne III is very unlikely, 
and these two lines are probably a consequence of the sky subtraction. Indeed, given 
that slit positions differ between observations, the slit during Obs 3 may be aligned with 
non-homogeneous structures, leading to background variations between the two nods. 

Finally, a total of three blue-shifted H I emission lines are present in the nod 1 and 
nod 2 spectra of Cygnus X-1 during Obs. 2 but not during Obs. 1 nor Obs. 3. Their 
velocities vary between —1600 to —2200 km s“\ which hints at an origin from the stellar 
wind of HD 226868. This is likely due to the orbital configuration of the system during each 
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observation. Indeed, Hanke et al. (2009) showed that one should expect the presence of 
blue- or red-shifted X-ray absorption lines from the stellar wind, depending on the orbital 
phase of the BH (Fig. 14 in their paper). Close to phase 0.5, these lines are almost “frozen” 
while they are maximally blue-shifted around phases 0.2— 0.3. During Obs. 1, Obs. 2, 
and Obs 3, the BH phases are 0.43, 0.27, and 0.58, respectively. Our results are therefore 
consistent with the prediction given in Hanke et al. (2009). 


6. Discussion 

The main results of our multi- wavelength study of Cygnus X-1 are; (1) the compact 
jets contribute significantly to the MIR emission, with a spectral break « 2.9 x 10^® Hz 
(Obs. 1), (2) when the jet carries more power, the cut-off frequency shifts as expected 
to higher frequencies and becomes undetectable in the MIR (Obs, 2), (3) the strength 
of the bremsstrahlung from the stellar winds is likely anti-correlated to the soft X-ray 
emission, which confirms the anti-correlation between the mass-loss and accretion rates (all 
observations) . 


6.1. Detection of the spectral break during Obs. 1 

This, to our knowledge, is the first measurement of the spectral break of a compact 
jet based on spectral fitting instead of photometry, and only the second detection of such 
a break in a BHXB after GX 339—4 (Corbel & Fender 2002). We can use the theoretical 
expression of the expected cut-off frequency, to see if our value is consistent with that of 
GX 339—4, scaled wdth respect to the BH mass, distance, and accretion rate. In the HS, 
if Fiy is the flux density at a given frequency y, M the accretion rate, and Mx the BH 
mass, we have for a radiatively inefficient jet, from Falcke & Biermann (1995) and Heinz & 



Sunyaev (2003): 
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M a 

^ 1 ^ <"> 

where is the cut-off frequency. 

For GX 339—4, we consider a distance of 8 kpc, a BH mass of 10 Mq, and a flux 
density at 15 GHz of 15 mJy, respectively (Zdziarski et al. 2004; Fender et al. 1997; Corbel 
et al. 2003). For Cygnus X-1 the numbers are 2 kpc, 11 Mq (Caballero-Nieves et al. 
2009), and 15 mJy, respectively. Using Eq 4, we find i^qx ~ 4.06 x i^cYG- where i/qx 
and i^cYG are the GX 339—4 and Cygnus X-1 cut-off frequencies, respectively. For our 
measured t'cYG ~ 2.9 x 10^^ Hz, this relation gives vqx ~ 1-15 x 10^^ Hz, which is in 
good agreement with the values derived in Corbel & Fender (2002). Therefore, despite the 
intrinsic uncertainties associated with the MIR spectral fitting, we believe that we see the 
cut-off frequency in the Cygnus X-1 MIR spectrum during Obs. 1. 


6.2. Implications for the jet geometry 

We can use the presence of the synchrotron break to constrain the jet geometry. As 
pointed out in Heinz (2006), the radio photosphere is located at a de- projected distance of 
about 5 x 10^^ cm from the black hole, based on the fact that about 50% of the flux are 
unresolved at an angular resolution of 3 mas in the VLBA observations reported by Stirhng 
et al. (2001). 

For a Blandford-Koenigl jet (Blandford & Koenigl 1979), the location of the 
photosphere along the jet, r=i is inversely proportional to the frequency, r=i,u oc. 

Given the marginally higher radio flux in our observation compared to the Stirling et 
al. (2001) observations, the radio photosphere should be located at distance of about 
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6 X 10^'^ cm. At a frequency of « 2.9 x 10^^ Hz, the IR photosphere would be located at 
a distance of 2 .^= 12 . 9 x 1013 hz 1-7 x lO’^^cm. Since the break from optically thick to thin 
occurs when the photosphere reaches the base of the jet, our result would suggest that the 
base of the jet has a scale of about half a light second. For a BH of mass M 10 Mq, this 
corresponds to about 10,000 gravitational radii, which is fairly large compared to size scales 
typically assumed for the base of the jet, usually assumed to occur on scales of a few tens 
to hundreds of gravitational radii. For example, the jet in M87 is resolved down to scales of 
less than 100 gravitational radii. 

A simple Blandford-Koenigl model is a severe oversimplification over the almost 
four orders of magnitude in frequency and thus photospheric scale between IR and radio. 
Deviations are almost certainly to be expected. However, any change would still have to 
obey the observed flat SED. As argued in Heinz (2006), simple geometric changes to the 
Blandford-Koenigl model do not alter the behavior of the jet significantly. For example, 
allowing the jet to be continually collimated (with 7?jet oc 2 ^) instead of being conical 
(C = 1) and requiring energy conservation such that p <x results in a slight altered 
dependence of 2r=i,i/ cx and a change of the spectral index from completely fiat to 

a = 1 — i (Heinz 2006). Bringing the base of the jet inwards (toward the BH) would require 
a (^ > 1, implying a steeper spectrum, while microquasars tend to have flat to slightly 
inverted spectra. A simple geometric explanation for the large implied scale for the base of 
the jet is thus not plausible. 

It is clear that the jet must be accelerated between the base of the jet and the radio 
photosphere, changing the Doppler parameter S along the jet. Since the optical depth to 
synchrotron-self absorption is proportional to P, a significant increase in 6 from the base of 
the jet, where the IR emission originates, to the radio photosphere, would boost the radio 
optical depth, effectively moving the photosphere outward by a factor of However, 
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given that the jet in Cygnus X-1 is likely not ultra-relativistic (Gleissner et al. 2004), it is 
unlikely that this could bring the base of the jet inward by more than a factor of a few. 

Finally, it is possible that the estimated resolved flux fraction of the radio of ~ 50% 
was overestimated in Heinz (2006). For example, if instead of 50%, only 25% of the flux 
were actually resolved in that observation, it would decrease Zr=o by an order of magnitude, 
again moving the size of the IR photosphere and thus the base of the jet towards scales 
more typically associated with jet acceleration. However, very long baseline observations 
tend to underestimate the resolved flux, rather than the unresolved flux. 

If the size scale over which the jet forms is really this large, it w'ould imply that large 
scale outflows (i.e., winds) from the disk are fundamentally involved in jet formation. 
Indeed, a recent study of GRS 1915+105 have shown that the accretion disk wind and the 
jets were competing for the same matter supply (Neilsen & Lee 2009) . Moreover, it has long 
been speculated that winds support jet collimation by providing lateral pressure support. 
This would make the jet in Cygnus X-1 qualitatively different from AGN jets which have 
been imaged down to scales significantly smaller than this in terms of gravitational radii. It 
would also imply that the jet originates over a region much larger than the X-ray emitting 
corona that is typically assumed to be the base of the jet. This underscores the need for 
further confirmation of the nature of the outflow responsible for the observed radio emission 
and, by extension, the large scale diffuse radio and emission line nebula around Cygnus X-1. 


6.3. Origins of the radio and X-ray/ -ray emissions during Obs. 1 

Fig 10 displays the contributions of the stellar photosphere (blue), the bremsstrahlung 
from the expanding stellar wind (magenta), and the compact jet extrapolated to the 
high-energy (green), from our fit of the Cygnus X-l’s SED in the case 2, scenario (a). 



21 - 


The simultaneous Ryle, Spitzer, and RXTE data are also superimposed, as well as the 
non-contemporaneous INTEGRALjl^MS data that were recently published (Laurent et al. 
2011). The authors detected for the first time polarized -ray emission beyond 400 keV, 
and they argue that it is due to compact jets. Our result is fully consistent with this 
scheme, as the optically thin part of the jet perfectly describes the INTEGRAL /IQl'Si data, 
providing strong evidence that we actually detect the spectral break around 2.9 x 10^^ Hz. 
Nevertheless, it is clear that the RXTE spectrum of Cygnus X-1 cannot be described by the 
compact jet alone, and one or several extra components are needed. The spectral similarity 
between the spectral indices of the IR and the X-ray spectra might suggest that one of 
these components is synchrotron-self-compton (SSC) emission in the base of the compact 
jet itself (Markoff & Nowak 2004; Markoff et al. 2005). However, this would require a 
large compactness of that region, in direct contradiction of our discussion regarding the 
implied large size of the IR photosphere and thus the base of the jet above. The more 
natural interpretation would be Comptonization of disk photons by the corona, independent 
of the optically thin synchrotron spectrum from the jet. Moreover, whether or not this 
result implies a disk origin of the X-ray emission in other sources like GX 339—4, whose 
NIR spectra are more consistent with a synchrotron origin of the X-ray, remains an open 
question. 


6.4. A clumpy stellar wind? 

Bremsstrahlung is found to be variable and to contribute significantly to the MIR 
emission of Cygnus X-1, more than the compact jet itself. The fact that its intensity is 
similar in the HS (Obs. 1 and Obs. 3), and lower in the failed-transition state (Obs. 2) 
is consistent with an anti-correlation of the mass-loss and accretion rates, as it was first 
shown by Gies et al. (2003), who derived a mass-loss rate of about 2.57 x 10“® Mq yr"^ 



and 2.00 x 10"® Af© yr"^ for the HS and SS, respectively. We can also derive from the 
theoretical expression which links it to our fit-derived bremsstrahlung flux density Sy at 
15 GHz, including clumping effects (Lamers &: Waters 1984; Scuderi et al. 1998): 


Mw = 2.26 10"^ X 


\/Tx 


( 10 GHz) 


- 0.6 


10^ K 


- 0.1 


D, 


1 kpc 


IJ-Voo 


V 100 km s 


( 5 ) 


where /oc is the clump filling factor of the wind, Te is the electron temperature far from 
the photosphere, which is well approximated as O.oT*. /r is the mean atomic weight of the 
expanding plasma per electron, usually taken at about 1.3 for O supergiant (Lamers & 
Leitherer 1993), and Uoo is the terminal velocity of the wind (« 1600 km s"^ Gies et al. 
2008). When clumping is not included (i.e. foo = 1), we derive mass- loss rates of about 
8.79 X 10"® Mq yr"^ and 5.95 x 10"® M© yr“^ for Obs. 3 (HS) and Obs. 2 (IS), respectively. 
Neglecting clumping clearly leads to overestimate the mass-loss rate by a factor of 3—4, as 
shown by Lamers & Waters (1984). This is therefore a hint that the Gygnus X-1 wind is 
clumped, and we need foo ~ 0.09 — 0.10 in Eq. 5 to match the mass-loss rates given in Gies 
et al. (2008). This value is consistent with that derived from the study of the X-ray dips 
(Hanke et al., in prep). 


7. Conclusion 

We have investigated, through SED fitting, the origin of the MIR emission of 
Gygnus X-1 when the source was in the HS (Obs. 1), the IS (Obs. 2), and in a peculiar 
compact jet- free HS (Obs. 3). We interpret the MIR continuum during Obs. 3 as due to the 
stellar’s photosphere and clumpy expanding wind (/oo 0.09 — 0.10). In contrast, when 
jets are present, they clearly contribute to the MIR emission. During Obs. 1, we detect the 
spectral break (i^b ~ 2.9 x lO’^® Hz), and we argue that this value is consistent with the 
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base of the jet being very large, likely involving accretion disk wind in the formation and 
collimation of the jet. We show that if the jet alone can account for the 7 -ray emission 
of the source, additional processes are needed to explain the 3 — 200 keV spectrum. A 
compact jet is also present during Obs. 2 but we do not detect the cut-off frequency as it 
has likely shifted towards higher frequencies. Nonetheless, the bremsstrahlung emission is 
found lower than that during Obs. 3, which is a confirmation of the anti-correlation between 
the mass-loss and accretion rates. 

Our results highlight the importance of the MIR and, by extension, the far-infrared 
in the study of microquasars, as all their components may have a spectral signature in 
these domains. The existence of powerful infrared space-based observatories, Herschel 
today and the James Webb Space Telescope (JWST) in the near future, combined to 
ground-based telescope such as the Atacama Millimeter/sub-millimeter Array (ALMA), 
are therefore a unique opportunity to make strong progresses in our understanding of their 
behavior and properties. We definitively recommend further studies of Cygnus X-1 and 
other microquasars with these facilities. 
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Table 1. Summary of all the observations of Cygnus X-1 we made use of in this study 


Obs. # 

IRS 

RXTE 

Ryle 

Ryle fluxes (mJy) 

4* 

1 

53513.0399839-53513.0544390 

53513.0433283-5351 3.2142542 

53513.0398560-53513.0546265 

5-28 

0.426 

2 

53528.9592175-53528.9736728 

53528.9698098-53529.2433283 

53528.9593506-53528.9737549 

4-42 

0.269 

3 

53553.112987) --53553.1274610 

53552. 9079579-53553. 198883S 

53553.1 131592-53553. 1275635 

0-48 

0.583 


Note. - We give the observation number, the IRS day of observation (in MJD), the day of simultaneous coverage with RXTE and the 
Ryle telescope, the Ryle flux level at 15 GHz (in mJy), as well as the orbital phase calculated from the ephemeris given in Brocksopp et al. 
(1999b)- 


OO 



Flux densit 





Table 2. Best parameters from the fits of Obs. 1, Obs. 2, and Obs. 3 RXTE /PCA+UEXTE spectra of Cygnus X-1 


Obs. # 

A^h 


infold 

Ti 

^brciik 

T2 

■^iron 

^iron 

F-i-20 

f’ 20-200 

X^/dof 


(:10“* cm”*) 

(keV) 

(keV) 


(keV) 


(keV) 

(kcV) 

(10~*nrg cm”* s”‘) 

(10""®erg cm“^ s"+) 


1 


26 1 5+2.09 

175.97l®;^J 

1.83l«;Sl 

10.24l«;i 

1 5^-fO.OI 

6 qo+0-06 

0-3K_o.o6 

0 60^ ® 
^•’^^-0.08 

].19 

2.41 

1.02/303 

2 

0.53 (frozen) 

26.871^2;;,^ 

151.67l^;i 

2.08i«:g; 

10.33l«;l| 

1 70+0.01 

^•^^-0.03 

O'Salo.ol 

1.55 

2.21 

1.08/304 

3 

1 04’^ 

25.281?,;“ 

212.18l«;“ 

1 7C + 0.01 

lOlllJl 

1 cq+0.01 

l.OO_o 01 

6 41 

'^•^^-0.04 

0 3Q+0 07 

0.71 

1.64 

1.09/303 


Note. — The bejst-fit model is phabs(highecutxbknpower+gaussian), and the errorbars are given at the 90% confidence level. 


‘^'Fhe fluxed are imabsorbed and exclude that of the iron feature 
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Fig. 4. — Hardness-Intensity Diagram of Cygnus X-1 corresponding to our three observations. The 
X-axis is the 1.2—12.0 kev ASM flux in Crab (1 Crab = 75.5 counts/s), while the Y-axis is the ASM 
HR2=C/A hardness ratio. Clearly, Obs. 3 is the faintest and hardest observation, while Obs. 2 is 
the brightest and softest. 









Table 3. Best parameters from the fits of the Obs. 3 MIR spectrum of Cygnus X-1 


Obs. # 

n 

Rt / R* Off 

Na 

X^/dof 


(K) 

(Hg/kpc) 

(inJy at 15 GHz) 


3 

28000 (frozen) 


A 't^-TO.O? 

1.00/302 


Note. ” ~ The model used is a combination of a black body and a power law. 
The errorbars are given at the 90% confidence level. 



Table 4. Best parameters from the fits of Obs. 1 and Obs. 2 radio/MIR SElDs of Cygiius X-1: case 1 


Obs. # 

% 

R*/D, 

Off 

Na 



xVdof 


(K) 

(HQ/kpc) 


(mJy at 15 GHz) 


(mJy at 15 GHz) 


1 

2 

28000 (frozen) 

9.94 (frozen) 

0.55 (frozen) 

0.66lg:“ 

0.44l«:“ 

o.ori°:“ 

0.06i“;“| 

15.601:®;" 

22.00lt:i 

0.96/301 

0.72/301 


Note. The model used is a combination of a black body and two power laws. The errorbars are given at 


the 90% confidence level. 



Table 5. Best parameters from the fits of Obs. 1 and Obs. 2 radio/MIR SEDs of Cygnus X-1: case 2, scenario (a) 


Obs. # 

(K) 

R,/D, 

(Hg/kpc) 

Off 

Ng 

(rnJy fit 1.5 GHz) 

Ojl 

Vc 

(THz) 

Qj2 

Ni 

(mJy at 15 GHz) 

X^/dof 

1 

2 

28000 (frozen) 

9.94 (frozen) 

0.55 (frozen) 

^ 7. + O.O 3 

rj 70+0.02 
^•'^- 0.04 

^^•'^-‘--0.02 
—0 01 

27.05lfjs 

14.80+:?,? 

-0.6 (frozen) 
-0.6 (frozen) 

15.40^3 + 

21 -30111 

0.81/300 

0.69/300 


Note. The model used is a combination of a black body, a power law, and a broken power law. The errorbars are given at the 90% 


confidence level. 


CaD 

■<I 



Table 6. Best parameters from the fits of Obs. 1 and Obs. 2 radio/MIR SEDs of Cygnus X-1: case 2, scenario (b) 


Obs. # 


R*/1X 

Off 

Nff 

Ojl 

1>C. 

Oj2 

fVj 

X^/dof 


(K) 

(Ii 0 /kpc) 


(sujy at 15 GHz) 


(T»i) 


(mJy at 15 GHz) 


1 

28000 (frozen) 

9.94 (frozen) 

0.55 (frozen) 

0.66 (frozen) 

o.o4l°:«l 

29.43l|g 


15.201°:“ 

0.85/300 

2 

= 

= 

= 

0.44 (frozen) 

0-06i«;«i 

42.6ltf:?| 

1 91+0.80 
0.36 

22.04l”;i 

0.69/300 


Note. The model used is a combination of a black body, a power law, and a broken power law. The errorbars are given at the 90% 


confidence level. 



Table 7. Best parameters from the fits of Obs. 1 and Obs. 2 radio/MIR SEDs of Cygnus X-1: case 2, scenario (c) 


Obs. # 

T, 

(K) 

R,/D, 

(Rlj/kpc) 

Off 

iVff 

(mJy at 1.5 GHz) 

aji 

Vc 

(■IHz) 

Qj2 

(mJy at 15 GHz) 

xVdof 

] 

2 

28000 (frozen) 

9.94 (frozen) 

0.55 (frozen) 

0.74 (frozen) 

nni+0-01 

A A1 +0.01 
0.01 

29 27+1|o 

1 T A T O. 44 
■ -^•-^'^-0.34 

_n r,r.4-0.04 

i5.3o:'Tg 

21 . 6011 :® 

0.78/300 

0.69/300 


Note. - I'he model used is a combination of a black body, a power law, and a broken power law. The errorbars are given at the 90% 


confidence level. 
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Fig. 6. — Fit of the MIR continuum of Cygnus X-1 during Obs. 3, with a stellar black body and 
an expanding wind bremsstrahlung. 
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Obs. 1 and Obs. 2 MIR/radio SEDs of Cygnus X-1: case 1 
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Fig. 8. — Obs. 1 and Obs. 2 MIR/radio SEDs of Cygnus X-1, case 2: (top) scenarios (a), (middle) 
scenario (b), (bottom) scenario (c). 











Table 8. List of all detected features for each LL 



position 


Nod 2 

A W FWHM Flux 

(;ini) (/'-m) W cm“^) 


16.20 -0.008i0.001 Q.24±0.01 0.96±0.16 

16.88 -0-019±0,004 0.37i0,02 1.22±0.23 

18-62 -0.022i0.001 0.2810.01 1.201,0.11 

19.10 -0.015±0.002 0-231.0.01 0.8610.06 


14.96 -0.00310.001 0.1110.01 0.3310.08 


19.10 -0.01010.006 0.0710.02 0.4710.30 


14.71 -0.01110.001 0.2610.01 1.3110.09 

16.20 -0.00810-003 0.3210.02 0.6410.16 

19.10 -0-04110,001 0.1810.02 1.741.0,06 


fuil-width at half-length (FWHM), and its flux, 
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Fig. 10. — Radio to X-ray SED of Cygnus X-1 during Obs. 1 (case 2a). The compact jet (green), 
bremsstrahlung (magenta), and stellar (blue) contributions are displayed. Simultaneous Ryle, IRS, 
and RXTE data are superimposed. We also add the data from the INTEGRAL/\^\S compton 
mode, obtained by combining all the IBIS observations of the source, in the HS, between 2003 and 
2010 (Laurent et al. 2011). 




